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Abstract.  Goldfarb  and  Hao  have  proposed  a  network  simplex  algorithm 
that  will  solve  a  maocimum  flow  problem  on  an  n-vertex,  m-arc  network  in 
at  most  nm  pivots  zmd  0{n^m)  time.  In  this  paper  we  describe  how  to 
implement  their  algorithm  to  run  in  O(nmlogn)  time  by  using  am  exten¬ 
sion  of  the  dynamic  tree  data  structure  of  Sleator  and  Tarjan.  This  bound 
is  less  tham  a  logairithmic  factor  larger  than  that  of  any  other  known  algo¬ 
rithm  for  the  problem. 
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1.  The  Maximum  Flow  Problem 
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Let  G  =  (F,  E)  be  an  imdirected  graph  with  vertex  set  V  of  size  n  and  edge  set  E  of 
size  m.  We  regard  each  edge  {u,ty}  as  consisting  of  two  oppositely-directed  arcs,  (v,w) 
and  {w,v).  For  any  vertex  v  we  denote  by  E(v)  the  set  of  vertices  w  such  that  {v,u>}  is 
an  edge.  We  assume  that  G  is  connected  and  that  n  >2.  Let  each  arc  (u,  u>)  of  G  have  a 
nonnegative  real-valued  capacity  u{v,w).  Finally,  let  s  and  t  be  two  distinguished  vertices 
of  G;  s  is  the  source  and  t  is  the  sink.  A  (feasible)  flow  on  G  is  a  real-valued  ftmetion  / 
on  the  arcs  satisfying  the  following  constraints: 
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For  every  axe  (v,u>),  f{v,w)=  —f{w,v)  (antisymmetry  constraints)  (11) 
For  every  arc  (v,tzj),  f{v,w)  <  u(r,ty)  (capacity  constraints)  (1.2) 

For  every  vertex  t;  ^  f{v,w)  =  0  (conservation  constraints)  (1-3) 

u)€£(ti) 

The  value  of  a  flow  /  is  value  (/)  =  fis,v).  The  maximum  flow  problem  is 

that  of  finding  a  flow  of  maiximum  value. 

To  date,  the  asymptotically  fastest  known  algorithms  are  those  of  Goldberg  and  Tarjan 
[8]  and  of  Ahuja,  Orlin,  and  Tarjan  [1].  The  former  nms  in  0(nmlog(n^/m))  time.  The 
latter  requires  integer  capacities;  it  runs  in  0(nm log(n(log  C7)^'^^/m  +  2))  if  no  capacity 
exceeds  U.  Both  of  these  algorithms  Eire  based  on  the  0(n^  )-time  algorithm  of  Goldberg 
[5].  Extensive  discussions  of  the  problem,  its  applications,  and  classical  algorithms  for  it 
can  be  found  in  [5],  [13],  [14],  [17]. 

The  above  statement  of  the  maximum  flow  problem  simplifies  notation  by  avoiding 
explicit  mention  of  “forward”  and  “backward”  residual  arcs.  It  is  completely  equivalent 
to  the  usual  formulation  on  directed  graphs.  The  case  where  two  oppositely-directed  arcs 
(vyw)  and  (u.-,u)  have  nonnegative  capacities  ti(r,u»)  and  u(u»,u)  and  zero  lower  bounds 
can  be  represented  by  an  undirected  edge  {y,  ty}  having  lower  boimd  —u{w,  v)  and  capacity 
u{v^w).  We  also  assume  for  simplicity  that  no  pair  of  vertices  in  G  is  connected  by  more 
than  one  edge,  but  allowing  G  to  be  a  multigraph  does  not  in  any  w’ay  affect  our  results. 

2.  The  Network  Simplex  Algorithms 

The  network  simplex  algorithm  is  a  specialization  of  the  revised  simplex  method  that 
uses  an  appropriate  data  structure  and  a  pivot  selection  rule  for  its  implementation.  It  is 
based  on  an  eEirly  observation  by  Fulkerson  and  Dantzig  [6]  and  Dantzig  [4]  that  any  basis 
matrix  of  a  vertex-edge  incidence  matrix  of  G  corresponds  to  a  rooted  spanning  tree  and 
can  be  permuted  to  eui  upper  triangular  matrix  with  a  ±1  diagoned.  (For  a  description  of 
the  method  see  e.g.  the  books  of  Chvatal  [2],  and  Kennington  and  Helgason  [12];  for  an 
implementation  see  Grigoriadis  [11].) 

We  state  the  network  simplex  algorithm  for  the  maximum  flow  problem  in  a  form 
suitable  for  our  implementation;  we  omit,  for  example,  the  return  arc  (<,s)  that  is  added 
in  the  standard  treatment.  Given  a  flow  /,  an  arc  (v,  u))  ha*;  rf$idual  capacity  Uf{v,ui)  = 
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u(r,  tr)  —  f(v,w).  Arc  {v,w)  is  saturated  if  Uy{v,w)  =  0  and  residual  if  Uf{v,w)  >  0.  An 
edge  {v,  u;}  is  saturated  if  either  (v,w)  or  (u>,  v)  is  saturated,  and  residual  otherwise.  A 
basic  flow  is  a  flow  /  such  that  the  set  of  residual  edges  forms  a  forest  (a  set  of  trees)  with 
s  and  i  in  different  trees.  Given  a  basic  flow  /,  a  basis  is  a  pair  of  trees  5,  Z  that  are 
subgraphs  of  G,  such  that  s  €  5,  <  G  Z,  and  every  vertex  and  every  residual  edge  is  in 
either  5  or  Z.  Given  a  basic  flow  /  and  a  basis  S,  Z,  an  edge  (or  arc)  is  a  tree  edge  (or 
tree  arc)  if  it  is  in  S  or  in  Z,  and  a  nontree  edge  (or  nontree  arc)  if  not.  A  basic  flow  /  is 
called  degenerate  if  there  is  a  saturated  tree  edge  and  nondegenerate  otherwise. 

The  network  simplex  algorithm  maintains  a  basic  flow  /  and  a  corresponding  basis 
5,  Z.  Starting  from  such  a  flow  /  and  basis  S,  Z,  the  algorithm  consists  of  repeating  the 
following  step  imtil  there  is  no  residual  arc  (v,u>)  with  u  G  5,  f  €  Z: 

Pivot.  Select  a  residual  arc  (v,iy)  wdth  v  £  S,w  £  Z .  .4dd  to  S  U  Z,  forming  a 

single  spanning  tree  T.  This  tree  contains  a  unique  simple  path  p  of  tree  arcs  from  s  to  t. 
Let  S  be  the  minimum  capacity  of  an  arc  on  p.  Add  S  to  the  flow  of  every  arc  on  p.  Delete 
from  T  some  edge  {x,t/}  such  that  (i,y)  is  a  saturated  arc  of  p.  This  produces  two  trees 
that  form  a  basis  for  the  new  basic  flow. 

Arc  {v,w)  is  called  the  entering  arc  of  the  pivot  and  (x,y)  the  leaving  arc;  the  pivot 
is  said  to  be  on  i,v,w).  It  is  possible  for  6  to  equal  zero  if  the  basic  flow  is  degenerate;  then 
the  pivot  is  said  to  be  degenerate.  A  degenerate  pivot  does  not  change  the  flow'  but  does 
change  the  basis.  A  nondegenerate  pivot  changes  the  flow,  increases  the  flow  value,  and 
may  or  may  not  change  the  basis. 

K  a  basic  flow  and  a  corresp  .i  ug  basis  are  not  available  initiailly,  they  can  be  com¬ 
puted  in  0{nm)  time  in  several  w  One  way  is  as  follows.  Let  /  =  0  and  compute  a 
spanning  tree  of  T  of  G.  Then,  select  a  nontree  residueJ  edge,  identify  the  unique  simple 
cycle  it  forms  in  T,  and  push  flow  around  this  cycle  so  that  at  least  one  of  its  edges  is 
saturated.  Repeat  this  step  until  there  are  no  nontree  residual  edges.  Finally,  push  enough 
flow  from  s  to  f  along  the  unique  (5,t)  path  in  T  so  that  at  least  one  additional  edge  is 
saturated.  Deleting  from  T  this  edge  yields  a  basic  flow  and  a  basis  5,  Z.  The  running 
time  for  this  computation  can  be  reduced  to  to  0(m  log  n)  by  using  the  dynamic  tree 
data  structure  [15],  [16],  [17],  but  this  does  not  improve  the  running  time  of  the  overall 
algorithm. 
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3.  A  Refinement  of  the  Algorithm  with  a  Polynomial  Number  of  Pivots 

The  algorithm  of  the  previous  section  need  not  terminate  unless  an  anti-cycling  rule, 
such  as  Cunningham’s  [3],  is  used  for  breaking  ties  in  selecting  the  leaving  arc.  For  integer 
data,  such  an  implementation  solves  the  maximum  flow  problem  in  at  most  n^U  pivots 
and  in  0{ri^mU)  time  using  a  simple  rooted  tree  data  structure  to  represent  the  basis. 
Goldfarb  and  Grigoriadis  [9]  proposed  a  mle  that  pivots  on  a  residual  arc  {v,xv)  with 
V  £  S,  XV  E  Z,  for  which  the  number  of  residual  arcs  in  the  paths  from  s  to  r  in  S  and  from 
uj  to  t  in  Z  is  minimum  over  all  nontree  residual  arcs  from  S  to  Z.  This  vauiant  works 
better  in  practice  than  others,  but  it  does  not  improve  the  pseudopolynomial  bound  on 
the  total  number  of  pivots. 

The  key  to  making  the  network  simplex  algorithm  nm  fast  is  to  choose  pivots  more 
carefully.  Goldfarb  and  Hao  [10]  proposed  a  pivot  rule  such  that  at  most  nm  pivots  occur. 
Explaining  their  rule  requires  a  few  extra  definitions.  We  call  an  arc  (v,  u))  pseudoreaidual 
if  it  is  residual  or  a  tree  arc^.  For  any  vertex  v,  we  define  the  label  d(v)  of  v  to  be  the 
minimum  number  of  pseudoresidual  arcs  on  a  path  of  pseudoresidual  arcs  from  s  to  v,  or 
infinity  if  there  is  no  such  path.  Every  vertex  label  remains  finite,  and  indeed  less  than  rt. 
until  after  the  l^^st  pivot.  Goldfarb  and  Hao’s  pivot  rule,  which  we  call  the  smallest  label 
rule,  is: 

Among  all  residual  arcs  (v,  i^)  with  v  €  5  and  u;  6  Z,  pivot  on  one  with  d{v)  minimum. 

Efficient  implementation  of  this  rule  requires  a  reformulation  of  it,  also  proposed  by 
Goldfarb  and  Hao:  Repeat  the  following  step  until  d(t)  =  oo: 

Choose  a  vertex  tv  ^  Z  with  d{w)  minimum.  Pivot  on  any  residual  arc  (v,  w)  with 
t;  €  S.  (Such  an  arc  will  have  d(v)  =  d(tv)  ~  1.) 

4.  Efficient  Implementation  of  the  Smallest  Label  Rule 

We  shall  describe  a  way  to  implement  the  smallest  label  rule  so  that  the  running  time 
of  the  resulting  network  simplex  algorithm  is  0(nm  log n).  This  improves  Goldfarb  and 
H21o’s  bound  of  O(n^m),  and  is  within  less  than  a  logarithmic  factor  of  the  bound  of  any 
other  known  algorithm. 

^  If  the  basic  flow  is  nondegenerate,  every  pseudoresidual  arc  is  also  a  residual  arc. 
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Oiir  implementation  consists  of  two  main  peals.  The  first  part,  described  in  this 
section,  is  a  way  to  maintain  vertex  labels  in  a  total  time  of  0(nm).  The  second  emd  more 
complicated  part,  explained  in  the  next  three  sections,  is  a  dynamic  tree  data  structure 
used  to  choose  pivots  and  to  maintain  the  basis.  The  amortized  time  ^  per  pivot  with  this 

data  structure  is  O(logn),  resulting  in  the  claimed  O(nmlogn)  overall  time  bound. 

• 

To  maintain  vertex  labels,  we  use  the  method  proposed  by  Goldberg  and  Tarjan  for 
maintaining  exact  distance  labels  in  their  maximum  flow  adgorithm  (see  [8],  Section  7). 
For  each  vertex  tv,  we  maintain  a  pointer  into  a  fixed  list  .A(ti))  of  the  arcs  (r,u)).  This 
pointer  indicates  a  pseudoresidual  arc  {v,w)  with  d{v)  =  d{w)  —  1.  That  is,  arc  {v,w)  is 
on  some  pseudoresidual  path  of  fewest  arcs  from  s  to  u).  We  call  (v,uj)  the  current  arc 
of  w.  For  each  vertex  w,  we  also  meiintain  a  list  L{w)  of  those  vei  tires  x  such  that  the 
current  arc  of  x  is  {w,x).  Initializing  this  information  at  the  beginning  of  the  ma-Kimum 
flow  computation  can  be  done  by  a  single  breadth-first  search  from  s,  taking  0{m)  time. 

Goldfenb  and  Hao  proved  that  vertex  labels  can  never  decrease,  only  stay  the  same 
or  increase,  as  the  .Jgorithm  proceeds.  Furthermore,  once  a  pseudoresidual  arc  {v,w) 
becomes  a  saturated  nontree  arc,  it  cannot  become  pseudoresidual  again  until  at  least  one 
of  d{v)  and  d{w)  increases. 

We  need  to  update  vertex  labels  after  each  pivot;  the  leaving  arc  (i,y)  may  no  longer 
be  pseudoresidual.  If  (r,y)  is  indeed  no  longer  pseudoresiduEil,  and  if  in  addition  (i,y)  is 
the  current  arc  of  y,  we  delete  y  from  L{x)  and  initialize  a  set  J?={y}  of  vertices  to  be 
relabeled.  Then  we  repeat  the  following  step  imtil  R  is  empty: 

Relabel.  Select  a  vertex  w  £  R  and  delete  it  from  R.  Let  {v,tv)  be  the  current  arc  of  w. 
(Since  tv  W2is  on  R,  either  (v,w)  is  no  longer  pseudoresidual  or  d{v)  >  d(w).)  Seem  the  arcs 
after  {v,w)  on  A{tv)  imtil  finding  one,  say  {x,w),  such  that  {x,w)  is  pseudoresidual  and 
d{x)  =  d{tv)  —  1,  or  reaching  the  end  of  A{tv).  In  the  former  case,  make  {x,w)  the  current 
arc  of  w  and  add  w  to  L(x);  the  relabeling  is  complete.  In  the  latter  case,  scan  all  of  A{w) 
to  find  the  first  pseudoresidud  arc  (y,  w)  on  A{w)  with  d(y)  minimum.  Make  (y,ie)  the 
current  arc  of  w,  add  w  to  L(y),  set  d{w)  =  d(y)  -f  1,  add  all  vertices  on  L{w)  to  R,  and 
set  L{rv)  =  0,  completing  the  relabeling.  If  there  is  no  such  arc  (y,tv),  then  d{iv)  =  oo; 

^  By  amortized  time  we  mean  the  time  per  operation  averaged  over  a  worst-case  sequence 
of  operations.  See  [18]. 
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the  maximum  flow  computation  is  complete. 

It  is  straightforward  to  verify  by  induction  the  correctness  of  this  method  of  maintain¬ 
ing  vertex  labels  and  current  arcs.  Each  arc  list  >l(iy)  for  w  ^  s  is  scanned  at  most  2n  —  2 
times,  twice  for  each  possible  value  of  d(w)  (from  1  to  n  —  1).  The  total  time  needed  to 
maintain  vertex  labels  is  thus  0(nm). 

5.  The  Use  of  Dynamic  Trees 

To  choose  pivots  and  maintain  the  basis,  we  use  an  extension  of  the  dynamic  tree 
data  structure  of  Sleator  and  Tarjan  [15],  [16],  [17].  This  data  structure  will  represent  a 
collection  of  vertex  disjoint  rooted  trees,  each  vertex  of  which  has  Jin  integer  label,  and 
each  edge  of  which  has  two  associated  reed  values,  g(v,w)  and  g(w,v).  We  denote 

by  pareni(v)  the  parent  of  vertex  v  in  its  dynamic  tree;  if  v  is  a  tree  root,  pareni(v)=nuU. 
We  adopt  the  convention  that  every  tree  vertex  is  both  an  ancestor  and  a  descendeint  of 
itself.  The  data  structure  supports  the  following  ten  operations  on  dynamic  trees.  Each 
operation  takes  O(log/c)  amortized  time,  where  k  is  the  toted  number  of  tree  vertices. 

make-iree{v)\  Make  vertex  t;  into  a  one-vertex  dynamic  tree.  Vertex  v  must  be  in  no  other 
tree. 

fiv.d-pareni(v):  Return  the  parent  of  vertex  v,  or  null  if  t;  is  a  tree  root. 

find-value(v):  Compute  eind  return  g(v,  parent(u));  if  y  is  a  tree  root,  return  infinity. 

find-min-value{v):  Find  and  return  an  ancestor  w  of  vertex  v  such  that  g(w^  parent(w)) 
is  minimum;  if  t;  is  a  tree  root,  return  v. 

find-min-label(v):  Find  and  return  a  descendant  w  of  v  that  has  minimum  label. 

change-label(v,  1):  Set  the  label  of  v  equal  to  /. 

change-value{v,S):  Add  real  number  6  to  g{w,  parent{w))  and  subtract  6  from 
g{parent(w),  tv)  for  every  nonroot  ancestor  w  of  v. 

link{v,  w,  a,  ^):  Combine  the  trees  contadning  v  and  w  by  making  tv  the  parent  of  v.  Define 
g{v,w)  =  a  and  g{tv,v)  =  Before  the  link  operation,  vertices  t;  and  tn  must  be  in 
different  trees,  with  t;  the  root  of  its  tree. 
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cu1(v):  Break  the  tree  containing  vertex  t;  in  two  bj  deleting  the  edge  joining  v  and  its 

parent.  Before  the  cut  operation,  vertex  v  must  be  a  nonroot. 

evert^v):  Reroot  the  tree  containing  vertex  v  by  making  r  the  root. 

To  implement  the  network  simplex  algorithm,  we  maintain  the  basis  S,  Z  as  a  pair  of 
dynamic  trees.  Tree  Z  is  permanently  rooted  at  <;  the  root  of  S  changes  as  the  algorithm 
proceeds.  Initialization  of  the  two  trees  requires  n  make-tree,  n  —  2  link,  and  n  change-label 
operations  at  the  beginning  of  the  algorithm.  Each  time  a  vertex  label,  as  computed  by 
the  method  in  Section  4,  changes,  we  perform  the  corresponding  change-label  operation. 

To  determine  w’hich  pivot  to  do  next  during  the  computation,  we  perform  find-min- 
label{t),  which  returns  a  vertex  in  Z,  say  w,  of  smallest  label.  We  pivot  on  the  current 
arc  {v,w)  of  w,  as  defined  in  Section  4.  To  actually  carry  out  the  pivc  we  first  perforn. 
everi{v),  to  root  S  at  v.  Then  we  perform  link{v ,  w ,  q ,  ^) ,  where  a  =  Uj{v,w)  2ind  (3  = 
u f{w,  v).  We  compute  the  leaving  arc  (i,  y)  of  the  pivot  by  letting  x  be  the  vertex  returned 
by  find-min-value{s)  and  then  letting  y  be  the  vertex  returned  by  find-pareni(x).  The 
amount  of  flow  to  be  moved  from  s  to  <  is  the  amount,  say  6,  returned  by  find-value{x). 
To  complete  the  pivot,  we  perform  change-value{s,  —6)  and  then  cut(x).  At  the  end  of 
the  maximum  flow  computation,  we  compute  the  flow  on  all  the  tree  arcs  by  using  n  —  2 
find  value  operations. 

With  this  implementation,  each  pivot  takes  0(1)  tree  operations.  The  amortized  time 
per  pivot  is  O(logn),  so  the  overall  running  time  of  the  network  simplex  algorithm  is 
0(nm  log  n),  as  desired. 

6.  Representation  of  Dynamic  Trees  by  Phantom  Trees 

It  remains  for  us  to  discuss  how  to  implement  dynamic  trees  so  that  the  amortized 
time  per  tree  operation  is  O(logn).  Obtaining  such  an  implementation  requires  extending 
the  Sleator-Tarjan  data  structure.  An  extension  designed  to  maintjiin  edge  values  and 
to  support  all  the  operations  except  find-min-label  and  change-label  appears  in  [19]  and 
can  be  used  without  modification  here.  The  novel  part  of  om  implementation  hes  in  the 
handling  of  vertex  labels;  whereas  the  original  dynamic  tree  data  structure  was  designed 
to  compute  combinations  of  values  over  tree  paths,  the  operation  find-min-label  requires 
combining  values  over  subtrees.  We  shall  describe  a  data  structure  that  supports  the 
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operations  make-tree,  find-parent,  find-min-label,  change-label,  link,  cut,  and  evert.  For  the 
other  operations,  we  can  either  use  a  sep2irate  data  structure  of  the  kind  described  in  [19], 
or  we  can  combine  the  two  structures  into  one.  This  can  be  done  by  adding  informaticr. 
representing  edge  values  to  the  data  structure  described  below.  (See  [19].) 

To  perform  find-min-label  operations  efficiently,  we  need  to  impose  a  constant  upper 
bound  on  the  valence  of  each  tree  vertex.  Thus  we  represent  each  dynamic  tree  D  by  a 
rooted  phantom  tree  P.  Tree  P  contains  all  the  vertices  of  D  and  possibly  some  additional 
dummy  vertices.  Each  vertex  in  a  phantom  tree,  henceforth  ceJled  a  p-vertex,  has  a  label  and 
a  color.  In  the  simulation  of  a  dynamic  tree  D  by  a  phantom  tree  P,  the  colors  are  vertices 
in  D.  Every  phantom  tree  has  maximum  valence  three.  The  following  operations  are 
allowed  on  phantom  trees:  make-tree,  find-parent,  find-min-label,  change-label,  link,  cut. 
emd  evert,  with  the  added  constraint  on  ItnA:  operations  that  link{v,w)  Ccinnot  be  performed 
unless  V  and  w  both  have  valence  at  most  two.  (The  third  and  fourth  paxcimeters  of  a  link 
operation  are  unnecessary,  since  edges  do  not  have  values  in  phantom  trees.)  Phantom 
trees  also  support  three  additional  operations; 

find-children{v):  Find  and  return  the  set  of  children  of  v. 

find-top{v):  Find  the  ancestor  of  v  closest  to  the  tree  root  that  has  the  same  color 

as  V. 

change- color{v,'))-.  Set  the  color  of  v  equal  to  7. 

The  precise  correspondence  between  dynamic  trees  and  phantom  trees  is  as  follows. 
In  a  ph2intom  tree  P  corresponding  to  a  dynamic  tree  D,  there  is  a  path  p{v)  of  vertices 
colored  v  corresponding  to  each  vertex  v  of  D.  One  of  the  vertices  of  p[v)  is  identified 
with  t;  and  has  the  same  label  as  v;  the  remaining  vertices  of  p{v)  arc  dummy  vert’cec, 
each  of  which  has  label  00  and  valence  exactly  three.  Each  edge  {u,u)}  of  D  corresponds 
to  an  edge  of  P  with  v'  colored  v  and  w*  colored  w.  That  is,  if  each  path  p{v)  in 

P  is  condensed  into  a  single  vertex  v  and  loops  (edges  of  the  form  {x,x})  are  deleted,  the 
resxilt  is  tree  D.  (See  Figure  1.) 

[Figure  1] 

We  simulate  each  of  the  dynamic  tree  operations  by  a  constant  number  of  phantom 
tree  operations,  as  follows; 
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make-iree{v): 

make-irte{v)\  change- coIot{v,v). 
find-parent{v): 

find-coloT{find-parent(find-top(v))). 

find-min-label{v): 

find-min-label{find-iop{v)). 

change-label{v): 

change-label{v). 

link(v,  w): 

Step  1.  Let  u=find-iop{v).  Peviorm  find- children^u).  If  u  has  two  or  fewer  children,  go 
to  Step  2.  Otherwise,  find  a  child  q  of  u  (if  any)  colored  v;  if  there  is  no  such  child,  let 
q  be  any  child  of  u.  Let  r  be  new  vertex  (not  in  any  phantom  tree).  Perform  make-tree 
(r);  change-color{r,v)\  change-label{r,  oc);  cut{q)\  lxnk{q,r)\  Itnkiu.r).  Replace  u  by  r 
and  go  to  Step  2. 

Step  2.  Perform  find- par ent{w)\  find-children{u').  If  u>  has  \'a]ence  two  or  less,  let 
X  =  w  and  go  to  Step  3.  Otherwise,  choose  a  child  y  (if  any)  colored  u ■;  if  there  is 
no  such  child,  let  y  be  any  child  of  u'.  Create  a  new  vertex  z  (not  in  emy  phantom 
tree).  Perform  make-tree{z)\  change-color{z,w)\  change-label{z,oc);  cui{y):  link{y.z): 
link{u\z).  Let  x  =  z;  go  to  Step  3. 

Step  S.  Perform  link{u,x). 

cut{v): 

Step  1.  Let  u=find-top{v)  and  x=find-parent{u).  Perform  cut{u).  if  u  =  u,  go  to 
Step  2.  Otherwise,  perform  find-chudTen\VL).  Let  q  and  r  be  the  children  of  u.  Perform 
cut(q)\  cut  (r);  link(q,r).  Destroy  dummy  vertex  u. 

Step  2.  If  z  =  tx;,  stop.  Otherwise,  find  the  two  vertices  y  and  z  adjacent  to  j  by 
performing  find-parent{x)  and  find-children(x).  If  one  of  y  and  z,  say  y,  is  the  parent 
of  X,  perform  cu<(z);  cut(x);  /infc(z,y).  Otherwise,  perform  cut{z)\  cut{y)\  link{z.y). 

In  either  case,  destroy  dummy  vertex  x  and  stop. 

I 
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eveTt{v): 

Each  dynamic  tree  operation  consists  of  0(1)  phantom  tree  operations  and  0(  1)  ad¬ 
ditional  work.  Since  each  dummy  vertex  in  a  phcintom  tree  has  valence  exactly  three  a 
dynamic  tree  containing  k  vertices  corresponds  to  a  phantom  tree  containing  at  most  3k'/2 
vertices. 

7.  Representation  of  Phantom  Trees  by  Virtual  Trees 

We  implement  phantom  trees  by  using  the  method  of  Sleator  and  Tarjan  [16],  modified 
only  as  necessary  to  deal  with  vertex  labels  and  colors.  We  assume  some  familiarity  with 
[IG];  w’e  shall  merely  sketch  the  details  of  the  implementation,  highlighting  the  changes 
needed  for  our  purpose.  (See  also  [17j,  Chapter  5.) 

We  represent  each  phantom  tree  P  by  a  rooted  virtual  tree  \\  w-hich  contains  the 
same  vertices  as  P  but  has  different  structure.  Each  vertex  of  V'  has  a  left  child  and  a  right 
child,  either  or  both  of  which  can  be  missing,  and  at  most  three  middle  children.  We  call 
an  edge  of  V  solid  if  it  joins  a  left  or  right  child  to  its  parent  and  dashed  otherwise.  Tree 
1’  consists  of  a  collection  of  binary  trees,  its  solid  subtrees,  connected  by  dashed  edges. 
The  parent  in  P  of  a  vertex  x  is  the  synunetric-order  successor  of  x  in  the  solid  subtree 
containing  x  in  V',  unless  x  is  last  in  its  solid  subtree,  in  which  case  its  parent  in  P  is  the 
parent  in  V  of  the  root  of  its  solid  subtree.  (See  Figure  2.)  That  is.  each  solid  subtree  in 
V’  corresponds  to  a  path  in  P,  with  symmetric  order  in  the  solid  subtree  corresponding  to 
the  order  along  the  path  from  deepest  to  shallow'est  vertex.  We  say  a  vertex  x  is  a  solid 
descendant  of  a  vertex  y  in  V,  and  y  is  a  solid  ancestor  of  y,  if  x  is  a  descendant  of  y  and 
the  path  from  x  to  y  consists  of  solid  edges. 

[Figure  2.) 

We  represent  the  structure  of  V'  by  storing  with  each  vertex  i  pointers  to  its  parent, 
its  left  and  right  children,  and  a  list  of  its  middle  children.  We  also  store  with  x  its  label 
and  color.  In  addition,  we  store  with  x  one  piece  of  cumulative  information,  min-label{x). 
which  is  the  minimum  label  of  any  descendant  of  x  in  V.  Finally,  we  store  with  x  a 
reversal  bit  rev(x),  used  to  handle  the  ct/ert  operation.  The  interpretation  of  reversal  bits 
is  as  follows.  Let  sum-rev{x)  be  the  mod-two  sum  of  the  reversal  bits  of  all  solid  ancestors 
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of  .  If  3um-rev{x)  is  1,  then  the  meanings  of  the  left  and  right  child  pointers  of  r  are 
reversed,  i.e.,  the  left  pointer  points  to  the  right  child  and  vice-versa. 

We  use  two  0(l)-time  restructuring  primitives  on  virtual  trees.  The  first  is  rotation. 
in  which  two  vertices  x  and  y  joined  by  a  solid  edge  are  interchanged  while  preserving 
symmetric  order.  (See  Figure  3.)  The  second  is  splicing,  in  which  the  left  child,  if  any.  of 
a  vertex  i  is  made  a  middle  child,  and  possibly  in  addition  some  middle  child  is  made  the 
left  child.  (See  Figure  4.)  A  splice  can  only  be  performed  if  i  is  the  root  of  a  solid  subtree. 
It  is  straightforward  to  verify  that  all  the  values  stored  at  each  vertex  can  be  updated  in 
0(1)  time  after  a  rotation  or  a  splice. 

[Figure  3] 

[Figure  4] 

The  main  restructuring  operation  on  virtual  trees  is  splaying.  A  splay  at  a  vertex  r 
consists  of  a  specific  sequence  of  rotations  eind  splices  along  the  path  from  i  to  the  tree 
root.  The  effect  of  the  splay  is  to  restructure  the  tree,  making  x  the  root.  The  actual  time 
required  for  a  splay  at  x  is  proportional  to  the  (original)  depth  of  x;  the  amortried  time  is 
O(logA-)  if  the  tree  containing  x  has  k  vertices.  See  [16]. 

We  can  perform  each  of  the  phantom  tree  operations  using  at  most  two  splay  opera¬ 
tions  and  0(1)  additional  restructuring  of  the  tree.  We  shall  describe  the  implementation 
of  three  of  the  operations;  implementation  of  the  oihcrs  is  similar.  (See  [16].)  To  perform 
tvtrt{v),  we  splay  at  v,  make  the  left  child  of  v  (if  any)  a  middle  child,  and  flip  the  bit 
Tev{v).  To  perform  find-min-label{v),  we  choose  a  vertex  x  of  minimum  mtn-label  among 
V  and  all  its  children  except  the  right  child.  We  search  down  through  descendants  of  x 
to  find  a  vertex  y  such  that  label{y)=min-label{x).  (This  search  is  guided  by  label  and 
min-label  values.)  Then  we  splay  at  y  and  return  y.  The  splay  at  y  pays  for  the  search  to 
reach  y.  We  perform  find-top{v)  as  follows.  First,  we  splay  at  v.  Then  we  let  t'  =  0  and 
t/Q  =  u.  We  repeat  the  following  step  until  n,  has  no  right  child  or  v,  differs  in  color  from 
V.  search  down  from  the  right  child  of  u,  through  left  children  until  reaching  a  vertex 
of  the  same  color  as  v  or  that  has  no  left  child;  replace  t  by  t  +  1.  Once  this  computation 
is  completed,  we  splay  at  v-  and  return  u-  if  it  has  the  same  color  as  u;  if  it  does  not,  we 
return  u,_i.  The  splay  at  v,  pays  for  all  the  searching. 

With  this  implementation,  the  amortized  time  per  phantom  tree  operation  is  O(log  n ). 
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This  implies  by  the  discmsion  in  Section  6  that  the  amortized  time  per  dynamic  tree 
operation  is  0(log  n).  By  the  discussion  in  Section  5,  this  implies  in  turn  that  the  amortized 
time  to  choose  a  pivot  and  implement  it  in  the  network  simplex  algorithm  is  0(log  n).  This 
gives  the  main  result  of  our  paper: 

Theorem  1.  The  Goldfaxb-Hao  version  of  the  primal  network  simplex  algorithm  for 
the  maximum  flow  problem  can  be  implemented  to  run  in  0(nm  log  n)  time. 
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Figure  3.  A  rotation  in  a  virtual  tree.  Triangles  denote  subtrees. 


